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We calculate two-photon exchange amplitude for the elastic electron-proton scattering in the 
framework of dispersion relations. The imaginary part of the amplitude is determined by unitarity. 
Since in the unitarity relation intermediate states are on shell, off-shell form factors are not needed 
for the calculation. The real part is then evaluated using analytical properties of the amplitude. The 
expression for the elastic contribution to the amplitude, obtained in our approach, differs from the 
results of traditional calculations with on-shell form factors. Nevertheless, numerically the difference 
is minor for Q^ up to 6 GeV^. 



I. INTRODUCTION 
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O ' The precision level of present-day electron-proton scattering experiments makes it necessary to take into account 

.^ , effects beyond Born approximation, such as two-photon exchange (TPE). TPE can be seen in various observables 
in wide kinematical range; in particular it influences proton radius measurements [l|, generates non-zero transverse 

^** , beam spin asymmetry [2|, and, the most important, TPE corrections play crucial role in reconciliation of different 

^^ ■ measurements of proton form factors (FFs) at high Q^ [3| . Clearly, such corrections are also required for analysis of 

I— i' data from upcoming measurements at higher Q^ 4]. 

,S^ • The TPE diagram (Fig. [T]) for clastic ep scattering differs from similar diagram in QED in two ways. First, the 

I ' proton is not a point-like object, thus there are some non-trivial FFs at 7p vertices. Second, the interaction of the 

' fj , proton with virtual photon may lead to excitation of inelastic intermediate states, such as irp, A resonance, and so 

^ ' on. 

Ch '. At present, the calculations exist for elastic intermediate state [3,0) and for a number of resonances Q^- Though 

'~~' • the evaluation of loop integral in these papers was almost perfect, the weak point of all such calculations is the starting 

, [ expression for the TPE diagram (here we consider the elastic contribution, but the contribution of resonances may 

h>. • be studied similarly). This expression results from the contraction of "Icptonic" and "hadronic" parts 
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where qi = k — k" , q2 = k" — k' , the leptonic part Lf^i, comes from QED 



_, k" + m 
^M- = ^ T^' 1772 Z:2 ^'^^ (2) 



but the hadronic part iJp,y is only guessed to be 



U\ ^-^ = ^'^'^^'^^^ ^^^ r.(9i)C/ + C/'r.(gi) |^±^ T,{q,)U (3) 




FIG. 1: TPE diagram 



^ There are also partonic model calculations, appropriate at large Q^ and energy Q. This approach will not be considered further. 



where p" — p + qi, p" = P + 92 and Tfj_{q) is the amphtude of proton interaction with electromagnetic field, written 
in the form 

r^('7) = l,F^{q^) - ±-F,{q')[j„q] (4) 

The above-described form of the TPE amplitude was used by many authors from Bodwin and Yennie in 1988 [9| to 
the latest papers [a, Q. The justification for such choice of i/^i/ is the following: first, it gives the expected result if 
the intermediate proton is on-shell, or more precisely, it has correct residues at p"^ = M^ and p"^ = M^, and second, 
this expression is gauge-invariant, i.e. 

qivHf,^ = q2t,Hf,^ = (5) 

However, (j31) is not the only expression with such properties. One can, for instance, add to F2 an arbitrary function 
which vanishes at p"^ — M^ , like 

F2{q^)^F2{q^) + {p"^-M^)f{q^) (6) 

We should emphasize that, if we are dealing with the elastic contribution only, than the choice of iJ^^ is somewhat 
a matter of convention, since any change of the elastic contribution may be compensated by appropriate redefinition 
of the inelastic one. Nevertheless, it is desirable to have clear, unambiguous, and easy-for-calculation definition for 
contribution of each intermediate state. 

The modification of 7p vertex at p"^ 7^ Af^, which of course can be more general than the example shown above, 
is usually referred to as introduction of proton off-shell FFs. The uncertainty of these FFs is believed to be the main 
source of theoretical uncertainty in TPE amplitudes |3[. On the other hand, such FFs are not directly measurable, 
just because off-shell proton cannot be a final state. Hence to take into account off-shell behaviour, one cannot 
rely on experimental data but instead must use some nucleon model. This is undesirable, since the result will be 
mo del- dependent. 

In the current paper we propose a consistent approach to calculation of TPE, in which the use of "off-shell" FFs is 
avoided. The approach is based on the dispersion relations. At first, the absorptive part of the amplitude is calculated 
using unitarity. Thus only "on-shell" FFs are needed to evaluate it. Then the whole amplitude is reconstructed by 
dispersion relations. Since this operation is linear, contributions from different intermediate states may be treated 
separately. 

II. THE AMPLITUDES 

We follow the notation of Refs.[i, [l^, [H. In particular, we define P = [p + p')/2, K = [k + k')/2 and t = q^, 
V — s ~ u ^ APK, where s, t and u are Mandelstam variables. The electron and proton masses are m and M, 
respectively. In the present section (but not in the whole paper) the electron mass is neglected. 

The general-case elastic ep scattering amplitudes is conveniently written as [l2| 

47ra - ( ~ 1 - P^-~\ 

M = -^u'l^u ■ U' ^^Fi - — [7^ q]F2 + J^KF.j U (7) 

In Ref. [101 the following set of amplitudes was introduced 



Ge = Fi- TF2 + iyFi/4.M^ 

Gm = Fi+F2+eiyF3/AM^ (8) 

G3 = i^F^/U'P 



which "diagonalizes" the cross-section 



da = ^^Y^ {e\GE\' + r|C?Mp + -^^^l^ap) (9) 

In the above equations, E is initial electron lab. energy, t = — i/4Af^ and e = [ly^ + t{4A-P — t)]/[i'^ — i(4M^ — <)]. 
Since the amplitude Gs vanishes in Born approximation and hence is 0(a), the last term in ([9]) is negligibly small and 
we have 

da^dao(e\GE\^ + T\GMn (10) 



similarly to Rosenbluth formula, except that Qe and Qm are £-dependent. 

However to make use of the dispersion relations, we need amplitudes, free from kinematical u and s singularities 
and zeros. Such amplitudes are easily constructed by consideration of annihilation channel. The helicity amplitudes 
of the process e~e~^ -^ pp are 



r++ = 47ra • 2i cos^ 9/2 {Jt{1+t)F:!, + F,„ + i^F^/UI^ 

r__ = Ana ■ 2i sin^ 9/2 (^/t{1 + t)F3 ~ F„, - i^F^/iM^ 

2M ( - ~ , , 

r+_ = r_+ = ATxa- —^ sin 9\F^ + vF-ijAM^ 

where Fe — Fi — TF2, Fm ~ Fi + F2, and 9 is i-channel scattering angle, 

cose = -i//v/-t(4A/2-i) 



(11) 
(12) 

(13) 
(14) 



The subscripts of the quantity T^^ indicate the signs of proton and antiproton helicities, respectively, while the 
electron and positron helicities are +1/2 and —1/2. Computing the scattering channel cross-section 



da 
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dt 647rAf2£;2 2 



{\T+ 



|T__p + 2|T+_p) 



(15) 



A+A+ii e 



(see e.g. 



we return to the formula ^. Each of the T^^ contains a kinematical factor of sin'^^^ ^' | cos 

Ref. [13|). The amplitudes free from kinematical singularities are obtained after removing these factors, i.e. 

v/r(l + t)F3 ± (Pra + uF^i/AM^-^ and F, + vF^/AM'^ 

or equivalently 

Gi = gE=F, + vF3/AM^, G2 = F„, + i^F3/4M\ G3 = F3 (16) 

The amplitudes Gn satisfy fixed-t dispersion relations 



1lGn{v) 



00 —l^th 



(17) 



and consequently, vanish at z^ — > 00. Under crossing i> — > ~i> two first amplitudes are odd and the last is even: 

GiM-^) = -Gi,2H, G3i-u)^G3iiy). (18) 



III. CALCULATION PROCEDURE 



A. Imaginary part 

The imaginary part of the scattering amplitude can be calculated via unitarity condition 

^if ^ -'- fi — Z^ ^ fn-^ in 



or graphically 



2Im 




h' h h' 

h, fi, „ „ h, 



z/Cn 



p 




(19) 



(20) 



where we have replaced T-matrix elements in the r.h.s. by their Born (one-photon exchange) approximations. Thus 
obtained is exactly the absorptive part of the TPE amplitude. 



Eg. ([20]) allows for natural and unambiguous classificatfon of different contributions to ImGn, according to inter- 
mediate hadronic states h. The term with h — proton will be called elastic contribution, the term with h — A(1232) 
will be the A resonance contribution and so on. Since the intermediate states appearing in the unitarity condition 
are real ("on-shell") particles, it is sufficient to know on-shell transition amplitudes of these states to calculate IraGn- 
Thus in particular the knowledge of proton "off-shell" FFs is not needed. 

The reconstruction of the Re G„ from Im Gn by dispersion integral is linear operation, therefore we may introduce a 
natural definition of elastic contribution to the whole amplitude as the quantity yielded by dispersion relation applied 
to the elastic part of Im G„ , and similarly for other contributions. 

B. Reconstruction of the real part 

From now on we consider the elastic contribution only. Such contribution to the imaginary part of invariant 
amplitudes G„ can be written in the form 

2 ^ _ 

ImG^r') = -|^ E / F^ih)F,it2)An,^,{ly,h,t2)9ik'^)Sik"^ - m^)9{p'i)Sip"^ - M^)d*k" (21) 

where Fi{t) = Fi{t)/{t — A^), An^ij is a polynomial in ti, ^2 and a rational function of ly (it may have poles in v; 
the explicit expression for An^ij is given in Appendix E]) . The 0- and S- functions ensure that intermediate particles 

are on-shell. The straightforward way of further calculation is to insert Im Gn into the dispersion integral p7|) and 
evaluate it. However this is not an easy task. For example, Ea. (j2ip for imaginary part is valid only in physical region 
1/ > \J—t{AM'^ — t), but the dispersion integral involves all v values above the threshold v — vth (corresponding to 
s — (M + m)'^), thus before it can be evaluated we must first find an analytical continuation of ((2T|) into the unphysical 
region. Though such analytical continuation is unique, it is hard to write it down in a compact form. So we will use 
an easier roundabout way. 

The amplitude G„ is an analytical function of ly with two branch cut discontinuities along the real axis: from — oo 
to — z^th and from z/th to +oo. As implied by Eq. ()17|) . it can be written as a sum of two parts, direct and crossed box 
amplitudes 

G„(j^) = Gn,hoA'^) + Gn,xhoA'^) (22) 

with each of them having only one discontinuity, box from z^th to -|-cx3 and crossed box from — oo to — i^th- Direct and 
crossed box amplitudes are related by 

Gn^box — iG„^xbox(~^) (23) 

where ± is chosen according to p^ . Thus to reconstruct G„ it is sufficient to find G„_box- 
To do this, we note that if we find any function with the following properties: 

1) it has no singularities except the branching point at s = (Af + m)^, 

2) its branch cut discontinuity is 2iImG„ , with ImG„ given by Eq. (f2T|) . 

3) it vanishes as s ^ cx), 

then such function necessarily coincides with the sought amplitude (otherwise their difference would be non-trivial 
bounded whole function, which is impossible). 
The analytical structure of FFs is such that 

oo _ 

in other words, the FFs in the Eq. pTjl are some linear combinations of a single poles j^ with a > 0. Using the 
decomposition (j24[) . we may obtain 

Q oc oo 

^ m,)F,{t2)A,,{u,h,h) ^Jdajdb ^^^ '^^'^";^^^_ ^) (25) 



and rewrite Eq. (l2T|) as 



ImG^:'^ = J da J dbc{iy,a,b) I j^—^^^—^^e{k'^)6{k"^ - m^)9{p'^)6{p"' - M^)d^k" (26) 

A2 A2 

Consider the function 

his, t- a,b)= J (^^ _ „)(^2 _ 5)(fc„2 „ ™2)(p//2 _ M^) ^^^^ 

It is well-known that this is an analytic function of s everywhere except the branch cut from s = [M + m)'^ to +00. 
Its discontinuity across the cut is 

-Mk'^)S{k"' - m^)9{p'^)S{p"' ~ M^)d'k" (28) 

[ti — a)[t2 — 0) 

which is exactly the innermost integral in (|26p . Thus if the coefficients c{i',a,b) were independent of v, the whole 
TPE amplitude would be obtained by substitution 

under the integral, yielding 

G^n-) - i^ Y.J m^)m2)An,,{-MM) ^^,„ _ J'l^l _ ^,^ (30) 



But actually quantities An^ij{v,ti,t2) and thus c{v,ti,t2) have poles at the boundary of the physical region v — 

±z/o = ±-\/— t(4M^ — t) (see explicit expressions in Appendix [A|) . Because of this function G„ , constructed by 
Ea. ([5n)l . satisfy conditions 2) and 3) but don't satisfy 1) since it has unphysical poles at v = ±vo. 

To remove these poles we may simply subtract the principal part of G„ Laurent series expansion about v = ±z^o 

N-l Af-1 

» - ^rn^) - E ^7:7^^;:^^ - E 7.7;7tB^ (31) 



^(cl) 



n.box^ J n w A. H (z. - ^o)^--^ ^^ r! (z. + ..o)^ 



where N is degree of the pole (actually 1 or 2) and 



9r±=g^{>^T>^or&,f{^ 



(32) 

^ = ±1^0 



Since the subtracted function is nicromorphic (has no branching points) and vanish at i/ = 00, the properties 2) and 
3) hold true and in addition, the obtained function Gn (v) is regular at 1/ — ztiyo- So the requirement 1) is also 
satisfied. Therefore G,j ^^^(z^) is the sought amplitude. 

In summary, the evaluation of the TPE amplitude proceeds as follows: 

1) construct the expression for the imaginary part in the form (|2ip . 

2) obtain the quantity (?„, Eg. ([50)1 . by substitution according to Eq. pS)) . 

3) subtract unphysical poles aX v — ±z^0: Eg. ((3T|) . 

4) perform (anti)synimetrization with respect to v, i.e. add crossed box amplitude. 

Due to decomposition (|25p the quantity G„ can be written as a linear combination of functions Ii{s,t;a,b) with 
different a and b. This is especially useful if FFs are parameterized as a discrete sum of a single poles (such an approach 
was used in Ref. p|). To perform the subtraction of unphysical poles one needs to know the value of function J4 
and its derivative at 1/ = zLi/q. They can be expressed via integrals similar to (|?7|) with k" — rn? or p"'^ — M^ or 
both dropped (such integrals were denoted /i, /2, h in Ref. @)- Some useful relations between them are given in 
Appendix [Bl With these relations, one may compare the expression for elastic part of TPE amplitude, obtained in 
the dispersion approach, with the "naive" result ^^. 
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FIG. 2: The amplitude change AGs 




FIG. 3: The TPE amplitude SGm obtained in old (dashed 
line) and new (solid line) approach. 



IV. RESULTS AND CONCLUSIONS 



After performing the above-described procedure, we have obtained the following results for the elastic contributions 
to the invariant amplitudes G„. The expressions for Gi and G2 remain the same as in the "naive" approach, Eqs.([T][3]) 



The expression for G3 is different: 



^ _ ^(naivc) 
(_rl — Lt]^ , 



G3 = G' 



G2 — G 



AG3(t) 



where 



AG3(t) = 



AttH 



m.)F2it2)(^^^^^^^^_ 2M. 



U t 



1 12 



fc"2 



d^k" 



The whole scattering amplitude may be written as 

47ra 



M = 7W(" 



a'ruU'iP^k - PKj^)U ■ AG3(t) 



g2jvf2 
Since the quantities that contribute to the cross-section up to the order 0{a) are 



Ge = Gi 



and 



Gm = G2 — 



4Af2 



(1 - £)G3 



(33) 



(34) 



(35) 



(36) 



(37) 



(see Eg. dTO)) ). with new expression ((34| for G3 TPE corrections to the cross-section will differ from those in "naive" 
approach, since 



m 



^(na 



=) _ ^r(l+T)Vl-e2AG3(t) 



(38) 



Moreover, the affected amplitude, Gm, is exactly the quantity which is responsible for the discrepancy between 
Rosenbluth and polarization transfer methods in the measurements of proton FFs [ll| . 

The numerical calculation, however, shows that the addition to Gm is very small (Figs. [5] and [31). Therefore most 
of the results obtained starting from "naive" expression for the amplitude will remain unchanged. In particular, we 
checked what the low-Q^ behaviour is the same as reported in Ref.[lO|, since the addition to Gm vanishes at Q^ -^ 0. 
Nevertheless, the proton off-shell form factors problem is overcome: they are not needed to calculate TPE amplitudes 
in our approach. 



APPENDIX A 

The coefficients An^ij may be computed in the following way. First, we write down the standard expression for the 
absorptive part of the amplitude (elastic contribution) 

Stt^ J qfq^ 

x0{k'^)S{k"^ - m'^)9{p'^)S{p"^ - M^)d^k" (Al) 

Then we decompose it into scalar invariant amplitudes according to Eq.([7]) (the equation for ImA^ contains Imi^„ 
instead of F„ in the r.h.s.). Since Im_P„ are scalars, they will depend on scalar combinations of p, k, p' , k', p" and k" . 
But due to " onshellness" of the intermediate particles we have p"^ = M^ and fc"^ = w? ^ and other scalar products 
can be expressed via g^, g|, v and t. The vertex functions F^ and Vy contain FFs, so the resulting expression will be 
quadratic in FFs. The ImGn are obtained as linear combinations of Imi^„, Eg. p^ . 

Below An a are written in a matrix notation, An = x"' ^"' ) • In these formulas tp = ti +^2 ^ t, t„i = ti — t2 

and vl = -i(4M2 - t). 

t{y-^{ (2t + tp)(4Af2 + ^-t) + 2Af% ) (^ 1 

■^ 16M2 1^^ + ^P + ^P ^2_ 2 [I 1 



ttpji^ - t){2MHp - i^{i^ - t)) ttit2\ /O 1 
16M2(z/2 _ 1,2) 4^2 J yi 2 

(t t(4M2 + j/-i) -ip(2M2 -t) ^iiM^ -t){2M^ + iy-t)\ f 2 1 

- (ly — tj <- + tp —-2 -jz h tt. 



2 " 2(i/2 - J.2) ' P (^2 _ ^2)2 J V 1 

"^2+ 2M2 + z/ - 1 1 / 



tt_^ (l + (^^^^-^ + V(--^) + ^^'^^^^ + 2(4M2 - t)(. - t)t 



16Af2 \^^ ^2_^2 --V-" WV-^ ^^^M^2 _ ^2)2 J ^ 2 



+ 



1/2 - vl "^"^ V 1 ; ' 4M2(zy2 _ J.2) ^^-^ \^ 2 

I _,., j tp{6M^ + i^-3t)- t{u - t) _ {SAP - t)v + (4Af 2 - t)(Af 2 -t)\ ( Q Q 
(y tjtp < /I/- 2 2\ i 



4(^2 _ ^2) P (Z.2_^2)2 J ^0 2 

^ ^ 2Af 2 + ,y-t fO \ , ^ ,^ , ^ , A/2(jy - t) / 1 \ , Af 2(,. _ t)tn,tp f -I 

-vtit2 5 2 — n 9 l+^pl^ + 'pJ' 



2/' P^ P^ ^2 _ j,2 I 1 / ' 7/2 _ j,2 MO 



+ (4A/2,,,,1M!±^ - 2A/2(. - 0^t,iM£!±f^%±i: 



(1.2 _ yl)2 J V 1 



2 1 



APPENDIX B 

It is convenient to use Breit frame, in which 



(A2) 



+ [tpii + ^p)4(^2_^2) - ^i*2| (^ 1 2 j + (^ - *)^™ ^ 2 ' 4(z.2 - ^2) ; ^ 1 

, (4M2-^)(:.-0^^ (2 l\ , ^(4Af2-t)(2A/2 + :.-t) ^^ (q o\ 

■^1^2 T n + 7772777?. ^2^ *i*2 no (A3) 



(A4) 



g-(0,0,0,V^), P-(iV4Af2-t,0,0,0), K ^ ^^=L== L, ^/^^^^,0,o] (Bl) 

At z/ = z^o a-ll three vectors have only time- and z- components. The components Pq and p" of the vector p" can be 
expressed via ii, ^2 and p"2 as 



The following identity holds 



s(ii +t2-t) 



(p"2 - M^){k"^ - m2) fc"2 - m2 



p"2 - RP 



(B3) 



(p"2-Af2)(fc"2 -m2) \ 4 2V4M2 - i 



Po 



v + AlVP -t^„ 



-Px 



For V = vq the r.h.s. vanishes. Multiplying the obtained equation by arbitrary function f{p") and integrating over 
d^p" we obtain the first sought relation 



fip")dY 



S{ti+t2~t) 



s~M^ 



s + NP - m^ 



(p"2-M2)(fc"2-m2) fc"2-TO2 



p" 



NP 







(B4) 



V = Vo 



To find the relations containing the derivative of /4, we divide Eq. (|B3p by 1^2 — ^2^ multiply by arbitrary function of 
the form f{p"'^,ti,t2) and integrate over df^p" , keeping in mind to put i/ — i/q afterwards. The r.h.s. will consist of 
three integrals, the last of which is 



f{p"\U,t2) 



I'rli^'l 



P'^d^p' 



^z/2 _ „2 J p"2_M^ s + p"^ + 2K^p'^ - 2{Po + Ko)p'^ 



(B5) 



(the long expression in the denominator is equal to fc"2 — ni^). Substituting p" — > — p" and averaging obtained and 
initial integrals, we obtain in the limit t/ — > j/q (which implies K^ -^ 0) 



0^23^ 



."^■i^" 



f{p"'MM) P'^d^P 



p"2 - NP fc"2 - m? 



f{p"\hM) P'i^d^p 



II2M„II 



y/AAP - t J p"^ - Af 2 (fc"2 - m2)2 

■'0 

After some algebra we obtain the second sought relation (the electron mass m was neglected in the numerators) 

s{ti+t2-t) s-A'P+m^ s + A'P -m'^} ti+t2-t 



(B6) 



'4^/' 



/(p"2,ti,t2)dV 



(p"2-M2)(fc"2-m2) A:"2-m2 



AP 



(B7) 



4j/t 



f{p"\hM)dV 



ti+t2-t 



ly^+t^ 



2vt 



A'P fc"2 - m? 



2vtit2 



2ttit2 



(fc"2-m2)(p"2-Af2) (fc"2_^2)2 
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